Stochastic heating of a molecular nanomagnet 
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We study the excitation dynamics of a single molecular nanomagnet by static and pulsed magnetic 
fields. Based on a stability analysis of the classical magnetization dynamics we identify analytically 
the fields parameters for which the energy is stochastically pumped into the system in which case 
the magnetization undergoes diffusively and irreversibly a large angle deflection. An approximate 
analytical expression for the diffusion constant in terms of the fields parameters is given and assessed 
by full numerical calculations. 
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I. INTRODUCTION 



In molecular nanomagnets (MNMs)i. such as in Mni2 
acetates the magnetic core of the molecule is surrounded 
by organic non-magnetic ligands that extinguish the in- 
ter MNM exchange interactions. Hence, much of the on- 
site physical properties are deducible by studying a sin- 
gle MNM, albeit the dipolar interaction is present and 
is essential for ordering phenomenal. Characteristic for 
MNMs is the relatively large effective spin (e.g., S = 
16 for Mni2) and the magnetic anisotropyi. MNMs ex- 
hibit a series of phenomenapi that are relevant for appli- 
cations in spintronics and quantum information^; most 
notably is the bistability behaviour, the resonance tun- 
neling of magnetization^ and the large spin relaxation 
time. The dynamical control and switching of the magne- 
tization via external fields is a key ingredient on the way 
to utilizing MNM for technological applications. In this 
context, the role of thermal and environmental effects 
have been considered^i^. For a single MNM^— "— at very 
low temperatures the main switching mechanism is the 
quantum tunneling of magnetization^. This is due to the 
large anisotropy barrier—. For an initially excited MNM 
and driven by external magnetic fields we have shown 
recentlyi^ that the phase space of the magnetization has 
a rich structure containing a separatrix of topologically 
different domains. Switching occurs at the separatrix 
as a consequence of a transition between these domains. 
However, the issue how to inject energy in this nonlinear 
system, i.e. how to realize the initially excited state near 
the separatrix, has not been addressed. Starting from 
the ground state, this question is not answered by res- 
onant fields as the system is non-linear, i.e. it changes 
its eigenfrequency as the oscillation amplitude variesi^. 
Formally the equations of motion for a single molecular 
magnet resembles the Landau-Lifshitz-(Gilbert) equation 
without the Gilbert damping. In the present case the 
negligible damping is an inherent system property and 
not a shortcoming of theory. This difference is insofar 
important as without dissipation a precessional switch- 
ing, e.g. as proposed inii, is not achievable. Hence, new 
switching schemes are needed for MNM that are different 



from those known for magnetic materials. One scheme 
proposed reeentlyi^ relies on a stochastic, diffusion-type 
switching. For this to work however, the system has to 
be excited to a desired state (near the separatrix). The 
question of how to achieve that is still open. In this paper 
we show that using appropriate polychromatic magnetic 
pulses we can achieve a stochastic heating of a molecular 
magnet as appropriate for the stochastic switching. We 
derive approximate analytical expressions for the field 
parameters that allow for the stochastic heating and test 
for our analytical predictions with full numerical calcu- 
lations. 



II. THEORETICAL MODEL 

We study a single molecular magnet (MM), e.g. Fes 
or Mni2 acetates and choose the z axis to be along the 
uniaxial anisotropy direction (easy axis). The MM is 
subjected to a constant magnetic field with an amplitude 
Hqi applied along the x-axis (hard axis) as well as to a 
series of magnetic pulses J-'{t) that are linearly polarized 
in the x direction. The Hamiltonian we write as^: 



H — Ho + Hi, 

Ho = -DS',+gfiBHoS^,Hi 



(1) 



giiBJ^{t)Sx 



Here D is the longitudinal anisotropy constant, Sx, Sy, 
and Sz are the spin operators projections along the x,y, 
and z directions, respectively, g is the Lande factor, and 
fiB is the Bohr magneton. Since the spin of the molecular 
nanomagnet is quite large a classical approximation is 
appropriate. Hence, it is advantageous to introduce the 
variables {Sz, ^p) via the transformation 



Sx = ^/l — SI cos(/3, Sy ~ \/\ — SI sin(/7 



and rewrite ([T]) in the compact fori 



,16 



H{t) ^ -~-Sl + VT~~Sf cos(p - VT~~Sf cos(^F(t) 



IDS 



,nt) 



gptsHa 



(2) 



Suppose that the apphed constant field is weak A > 1 
and at the initial moment of time the system resides near 
to the ground state Sz(t = 0) w ±1. How to pump ef- 
ficiently energy into the system such that we reach the 
excited states near the separatrix, where then a dynam- 
ically induced switching is realizable? To answer this 
question we make a transition to the action-angle vari- 
ables (/, ip) in which the Hamiltonian ^ reads 



H = Ho{I) + V{I,^)F{t), 



HO^LO{I)I,LO{I) = 



d/(I]) 



(3) 



/(S) = i /"5,(S,^)dv?, ^ = -H. 



The equations of motion (EOM) are 



I = -- 



dH 



dV{I,ip) 



dip 



F{t), 



dip 
dV{I,p) 

di 



Fit). 



(4) 



In absence of the time dependent perturbation, / is an 
integral of motion {ip{t) is fast variable, however). For 
an applied monochromatic field it is not possible to keep 
in resonance with w, for cj depends on / and hence it 
changes in time. A polychromatic field offers a wider 
range of frequencies that may match the dynamical fre- 
quency of the system. To be more concrete let us assume 
the applied field to be of the form 

F = eoT Y, Sr{t-nT) = eo^ Y. cosU^ni , 



-1/t 



(5) 

where t is the pulse duration, T > r is interval between 
pulses, and Sq is the pulse strength. 



dl 

duj{I) 



di^{I)dY. _ .duj{I) 



dY. 2(S 

Thus we infer 



3)^2 



■ E-l^ 



-fi = £a-L :tz _.. ^ — n — ^ > t- 



(S - 3)/v 3 



(8) 



vS-l; 



E 



jT3y); A'(j^^) are the complete elliptic integrals in 
the notation of Ref . [l^ . In the regime of chaotic motion, 
when Eq.(6) holds, a dynamical description becomes in- 
appropriate. The adequate language for the study of the 
magnetization dynamics in this case is an approach based 
for example on the Fokker-Planck equation. A Fokker- 
Planck equation for the distribution function of the ac- 
tion /(/,i) can be set up in a similar way as done in 
Rei.m-. 



at 2di ^' di ' ^' 



'T. 



X 



MI)). 



The relevant qua ntity is the averaged value of the action 
(/(i, S))/ = I{t). Using the relations 



dfil,t) 
dt 



di^U4Dii)'-i^di, 

2J dl ^ ' dl 



(9) 



n,.-i/.„M^.4 



dD{I) _ elT^ 



di 



2A 



uj{I)[l + I 



dD{I) 
dl 



dijj{I) 



fil,t)dl, 



(10) 
(11) 



III. STOCHASTIC HEATING 

Of particular interest for us is the situation of overlap- 
ping resonances which is realized wheni^ 



K' = eor/(S) 



dw(/) 

di{j:) 



>i, 



(6) 



in which case the dynamics turns chaotioi^,, i.e. the sys- 
tem jumps from one resonance to other in a random way. 
A key point here is the irreversibility of the dynamics 
that emerges due to nonlinearity and without any ther- 
mal effects nor external random forces. Hence, we expect 
a "stochastic heating" of a MM subjected to the pulses 
dl]) when the criterion (JH) is fulfilled. 

Assuming that ^^(i = 0) « ^1 and A > 1 we find 



/(E) 



^(I) 



S,X'Z.,p)dp 
dI(Y) 



Vx 

VA\/e 



VE- IE 



E-1 



Kilh) 



(7) 



we infer that 



2A 



i^{I)[ l + I 



dLj{I) 

dH 



(12) 



We note that / and w(/) are functions of E and make the 
approximation that 



/(E) = /(E), Q{I)^u{l{t)), 



(13) 



meaning that correlations between the random variables 
are neglected. Formally, we can systematically improve 
on this approximation by accounting for higher moments 
for the correlation functions of the random variables. Ne- 
glecting these correlations we find for / 



t=4i^--f^^/)fi+/^Y 



dEdt 2A"(^H' + ^^^J' 



dE 
Itt 






-'^) 



(14) 



(15) 
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FIG. 1. A demonstration of the stochastic deflection of the 
magnetization vector from the initial state Sz(0) = 0.9 upon 
applying a series of rectangular pulses. The results are ob- 
tained by solving numerically for the Hamiltonian equations. 
The time scale is set by values of the constant magnetic field 
t — >■ t/g^sHo [16]. The parameters are eo = 0.4, T — Ips, 
T = O.Olps, A — 10. Note, the criterion of the stochasticity, 
given by ((8]) is realized since K' = 1.75 > 1. 



we conclude that 



— - -siT- 



1 



1 



^From the asymptotical solution of S > 1 






m 



A 



■ exp 



-elTt 



(16) 



(17) 



we uncover a diffusive decay of S(i), meaning that the 
energy is increased diffusively due to the relation H = 
— S(i)), albeit eq. (|8]) must be obeyed. 



IV. NUMERICAL RESULTS 

The exact numerical results shown in Figs.(ITJ[5J[3]) evi- 
dence that the mechanism of stochastic heating is indeed 
present and is quite efficient. 

Other scenario for the magnetization control is to cm- 
ploy only the periodic scries of rectangular pulses ap- 
plied along the hard axis, i.e. to switch off the static 
field. Measuring the energy in units of D we write for 
the scaled Hamiltonian 



H{t) = H{t)/D 



5? 



+ V{s„^), 



where 



+ 00 



V{s,,^) = V„is,,ip)T Y, Sit-kT), 



/C— — CXD 



and 



Vo{sz,ip) = ea/I - slcosifi, 




tCps) 



FIG. 2. Illustration of the stochastic heating by exter- 
nal driving fields after having started from the initial state 
5z(0) = 0.9. The results are obtained as in Fig. [T]with the 
parameters eo ~ 0.4, T — Ips, t = O.Olps, A — 10. The 
criterion of stochasticity ([8]) is fulfilled (K' = 1.75 > 1) and 
a diffusive increase of the system energy is observed. This 
numerical result is qualitatively consistent with the analyt- 
ical prediction, given by Eq. (17). For a better numerical 
agreement one has to go beyond the approximation (13) and 
consider higher moments for the correlation functions of the 
random variables 



0.95 



0.85 



S(t) 0.8 




0.75 



0.65 



20 40 60 SO 100 120 140 160 ISO 200 
t(ps) 

FIG. 3. The time evolution of the magnetization vector when 
starting from the initial state Sz{0) = 0.9. The results are 
obtained by solving numerically for the Hamiltonian equa- 
tions with a series of rectangular pulses. The parameters are 
eo = 0.1, r = Ips, r = O.Olps, A = 10, K' = 0.43 < 1. 
The criterion of stochasticity is not fulfilled. The orientation 
of the magnetic vector is not changed significantly, and only 
fiuctuations around the ground state is observed. 



with 






The equations of motion 



ip = uj{sz) + e 



9Vb(s2,V3) 



+ CX3 



ds. 



T J2 Sit-kT), 



(18) 
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FIG. 4. S^{t) as a function of the number of applied kicks 
t = nT obtained from a full numerical integration of the re- 
currence relations ea. (|19p . The basic control parameter is 
the time interval between the pulses T which allows a tuning 
of the values of of stochasticity coefficient (because Ddtff ~ 
3(13 ) ^) ^^^ hence allows a realization of different types of 
dynamics. We consider the following initial conditions for 
JIl}: a) e = 0.1, &(0) = 0.99, v{0) = 0, T = O.Olps, and 
iCo = 10"^ b) e = 0.1, 5^(0) =0.99, v:'(0) = 0, T = O.lps, 
and A'o = 10"^ 



can be integrated exactly in this case by formulating 
them as recurrence relationsi^ using the evolution op- 
erator T that propagate the system from the time to to 
T, i.e. 

{JI,^) = T(s2,(/3),s7= Szih + T - 0), 



(p = ip{to + T - 0), Sj = Sz{to - 0), (/3 = ip{to - 0). 

Note, that T = TjiTg consists of two parts, one describ- 
ing the free rotations Tr and the other, Ts, the action of 
the applied pulses, i.e. 

fR{sz;(p) = {sz;(p + uj{sz)T). 
For Ts virc find upon integrating EOM 



Ts{sz,ip) = s^ - e ,tp + eT- 



dip 



Thus the recurrence relation applies 



dip 



sTyJl 



SJ Sin(y9, 



(19) 



Vn = iPn 



,T + erVl-s?sin^- 



eTsz 



V^ 



: COS (p. 



Depending on the chosen parameters, these relations ()19p 
may be be stable or unstable as signified by the cor- 
responding Lyapunov exponents. Here, we inspect the 
Jacobean matrix 



M 






and find for the eigenvalues 



K = eT 



Ai,2 = l + ^A-± J(l + ix)2-i, 



2 52^0 



d(p'^ 



zT'^\/\ — si cos tp = Ko\/l — si cos tp, 



where Kq = maxK ~ eT^. Chaos is expected if Ai > 1, 
i.e. for K > 0, meaning even weak pulses Kq = eT"^ > 
may lead to a diffusion. Lc., S^{t) = S'{t = nT) and 
the magnetization can be deflected diffusively and ir- 
reversibly if Kq exceeds a critical value K^, as demon- 
strated by the numerical calculations in Fig|3l Essential 
for this phenomena is the existence of two time scales, the 
slow variables S^{t) and the fast random phase (p{t). Us- 
ing the random phase approximation for the fast phases^S 
one infers the Fokker-Planck equation 

dt ^2ds, ^''' ds, ' 



where 



Dis.) 



Tie- 



2 +°° 



Y. m^\Vm{s,,)\^ n 



m— — cxD 



27r 



and Vm{sz) are the Fourier coefficients as deduced from 
the expansion 



+00 



yoisz,'fi)= Y Vjn{sz)exp{inif). 

n— — 00 

Explicitly we have 



Kf^v 



(e''^ + e""^). 



Therefore, 

Vi == F_i = - 
and the diffusion coefficient is 



Disz) = ^{\V-,\' + \Vin ^ ^i^-4) = "-^(^-4)- 



V^-^l 


2 


V^-sl 



Consequently, we write 






dt 



where 



D, 



'T 



diff 



is the coefficient of diffusion which is completely defined 
by the pulse parameters 

+1 

j^{Sz) = -2D{S,), (S,) = J Szf{Sz,t)dSz. (20) 



Thus, the mean value of the spin projection behaves as 



2-Ddiff* 



(5.(t)) = (5,(0))e-^^-// 



the analyticaUy estimated decay rate 1/Ddiff coincides 
with the numerieaUy deduced one (cf. FigH]). 



SUMMARY 



Since the diffusion coefficient is given by 



Dd, 



ff 



4 



D 



T 



{D is the magnetic anisotropy constant) Ddiff can be 
tuned by changing appropriately the external field pa- 
rameters, e.g. by varying the amplitude of the pulses Eq 
and/or the interval T between them. Depending on these 
parameters different types of the dynamics is realized. 
To test for this analytical prediction we performed full 
numerical calculations that are in good accord with the 
analytical results (see Fig|4]). This statement is based on 
the fact that for e = 0.1, T = 0.1, Dd^ff = e^T = 0.001 



The aim of this work is to point out the possibility of 
a stochastic energy pumping and magnetization defiec- 
tion in a single molecular magnet subjected to a static 
and a time-variable, polychromatic magnetic fields. The 
key point is that the parameters of the applied static 
and pulsed magnetic fields can be tuned such that the 
system is driven nearby a separatrix where the magneti- 
zation dynamics turns diffusive allowing thus for a mag- 
netization switching even in the absence of damping (that 
conventionally originates from coupling to other degrees 
of freedom). 
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